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Abstract 

A structural transformation in water upon compression was recently observed at 
the temperature T = 277 K in the vicinity of the pressure p ~ 2 000 Atm [R.M. 
Khusnutdinoff, A.V. Mokshin, J. Non-Cryst. Solids 357, 1677 (2011)]. It was found 
that the transformations are related with the principal structural changes within 
the first two coordination shells as well as the deformation of the hydrogen-bond 
network. In this work we study in details the influence of these structural trans- 
formations on the vibrational molecular dynamics of water by means of molecular 
dynamics simulations on the basis of the model Amoeba potential (T = 290 K, 
p = 1.0-7-10 000 Atm). The equation of state and the isothermal compressibility are 
found for the considered (p,T)-range. The vibrational density of states extracted for 
T H z-irequency range manifests the two distinct modes, where the high-frequency 
mode is independent on pressure whereas the low-frequency one has the strong, 
non-monotonic pressure-dependence and exhibits a step-like behavior at the pres- 
sure p ~ 2000 Atm. The extended analysis of the local structural and vibrational 
properties discovers that there is a strong correlation between the primary struc- 
tural and vibrational aspects of the liquid-liquid structural transformation related 
with the molecular rearrangement within the range of the second coordination shell. 
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1 Introduction 



The physico-chemical features and soluble properties of water are of the cen- 
tral importance for life on our planet. However, despite of a knowledge of 
some properties at the ambient conditions that corresponds only to narrow 
part of water phase diagram, the inherent structural and dynamical proper- 
ties are still far from to be completely understood, that motivates to per- 
form a large number of experimental and theoretical investigations of water 
[1,2,3,4,5,6,7,8,9,10,11,12,13]. Water is known as a 'system' with a rich variety 
of anomalous behavior related with the different phase states. For example, in 
a solid state, water has at least fifteen crystalline forms, where four of these can 
coexist with the liquid phase, and a variety of the amorphous phases known 
as amorphous ices and polyamorphism [2,14,15]. 

One of the interesting features of water is related with the first-order like phase 
transition from low-density liquid (LDL) state to high-density liquid (HDL) 
state was reported in Ref. [16]. The hypothesis of a "liquid-liquid phase tran- 
sition" in water was considered in a series of experimental, computational and 
theoretical studies (see, for example, Refs. [17,18,19] and references therein). 
As it was found, the difference between these liquid states is mainly provided 
by the structural properties. Namely, in the HDL, the local tetrahedrally coor- 
dinated hydrogen bond (HB) structure is not fully developed, whereas in the 
LDL, a more open, locally "ice- like" HB network is realized [1,3]. Both these 
liquid states are closely associated with the two amorphous phases observed ex- 
perimentally in water: the low-density amorphous (LDA) and the high-density 
amorphous (HDA) ices, although locus of this "transition" is still debatable 
[15,18,19,20]. Remarkably, that a sharp jump in density (volume) is smoothed 
with the increase of temperature (or the decrease of pressure) and disappears 
completely, while LDL/HDL structural transformations still exists. It is in- 
teresting to note that the amorphous phase of water is also characterized by 
the "transition" from HDA ice to very high-density amorphous (VHDA) ice 
[21,22,23,24,25]. A schematic phase diagram of water, that treats qualitatively 
some empirical results, is given in Fig. 1. Above the critical temperature of the 
predicted LDL/HDL transition, the jump in volume (density) is absent, but 
the local structure changes completely with pressure. This is clear detected 
by the radial distribution function (for example, for oxygen atoms), goo{ r ), 
where the first minimum disappears, while the second maximum and mini- 
mum changes by a minimum and maximum, respectively (see right panel in 
Fig. 1). It is significant to note that although a schematic phase diagram of 
water is known in general sense for an extensive (p,T)-range (see, for example, 
Fig. 1 in Ref. [25] and Fig. 4 in Ref. [26]), the exact boundary location between 
LDL and HDL is not well-defined as well as the smoothed LDL/HDL trans- 
formations range are still undetected. So, the problem associated with the 
structural transformation and its influence on microscopic dynamics of liquid 
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Figure 1. (Color online) Left panel: Schematic phase diagram of liquid-amorphous 
water (see also Ref. [27]). White range corresponds to a stable liquid water, gray 
region indicates the liquid-amorphous phase, which is characterized by the sepa- 
rated low-density and high-density states. Broken line is the conventional boundary 
between the liquid and the amorphous phases. Solid line is the 'isotherm' considered 
in the work, where a sharp jump in density is no observed. LDL/HDL-structural 
transformations arise at the temperature T € [260 K; 400 K] and the pressures 
p € [2000 Atm;4000 Atm] (see Fig. 5 in Ref. [28] and Fig. 7 of Ref. [29]). Right 
panel: Radial distribution functions goo( r ) typical for three ranges of the phase 
diagram (LDL, LDL/HDL-boundary, HDL); the correspondence is indicated by ar- 
rows. 



water (even at the room temperature) requires the additional studies [15]. The 
given work studies the structural features within first two coordiantion shells 
[30] as well as the vibrational properties of liquid water at the temperature 
T = 290 K and a wide range of pressure by means of molecular dynamics 
simulations on the basis of the Amoeba interaction potential. The correlation 
between structural features of the transformation and molecular vibrational 
dynamics of water is also considered. 

This paper organized as follows. The details of molecular dynamics simula- 
tions of liquid water are presented in section 2. The structural transforma- 
tions within the second coordination shell are explored by means of the local 
isothermal compressibility in section 3. The manifestation of the effect of 
the structural transformations on the vibrational properties of water and the 
possible correlation between these both are also discussed here. Finally, we 
conclude with a short summary. 
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2 Simulation Details 



Equilibrium classical molecular dynamics (MD) simulations of water at the 
constant temperature T = 290 K and the different pressures were performed 
on the basis of the AmoebaPI (Atomic Multipole Optimized Energetics for 
Biomolecular Applications) water model suggested recently by Ren and Pon- 
der [32,33]. The model uses a polarizable atomic multipole description of 
electrostatic interactions, where the polarization is treated due to the self- 
consistent induced atomic dipoles [32], and a modified version of the Thole's 
interaction model is used to damp the induction at a short range [34,35]. Our 
computations were performed for 4000 water molecules interacted within a 
cubic box with the periodic boundary conditions in all directions. The Ewald 
summation was used to handle the electrostatic interactions. Moreover, an 
atom-based switching window of the size 12 A was applied to cut off the van 
der Waals interactions. The equations of motion were integrated via a modified 
Beeman algorithm [36], 



f(t + At) = f(t) + $At + ^ ( Aa{t) - a(t - At) ] At 2 + C(Ar 4 ), 
$(t + At) = tf(t) + ^[ 5a(t + At) + 8a(t) - a(t - At) ] At + 0{At^ 

with the time step At = 1.0 fs. The isothermal-isobaric ensemble at the tem- 
perature T = 290 K and the pressure p = 1.0, 1000, 2000, 2500, 3000, 3750, 
5000, 7797 and 10 000 Atm was applied by means of the Berendsen ther- 
mostat and barostat [37]. According to this scheme, to enforce the constant 
temperature the system was weakly coupled to a heat bath with some tem- 
perature. The velocities were scaled at an each step by such a way, that the 
rate of the temperature change is 



dT = Tp-T 
dt t 



1 Recently, it was shown in Ref. [1], that the results of molecular dynamics simu- 
lations on the basis of Amoeba interaction potential for structural and dynamical 
properties of water are in a good agreement with the results on ab-initio molecular 
dynamics simulations as well as the experimental data on neutron diffraction [31] 
and inelastic X-ray scattering [6] . Note that the Amoeba model reveals for the isobar 
P = 1.0 Atm the maximum in density at the temperature T = 290 K instead of the 
well-known value 277 K [32]. Actually, one notices that the Amoeba potential does 
not reproduce the experimental equation of state correctly: for instance, the density 
1.21 gjcrr? at 290 K is reached at 5000 Atm in the present simulation, instead of 
8000 Atm in the experiment. 
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where r is the coupling parameter, which determines how tightly the bath and 
the system are coupled together, T is the temperature of the external heat 
bath. This method provides an exponential decay of the system temperature 
towards the desired value. The change in temperature between the successive 
time steps is 

AT=^(r -T(t) y j, (2) 

here 

1 N p 2 

T{t) = (3N-N c )k B § m, (3) 

the instantaneous value of the temperature; N c is the number of constraints 
and (3N — N c ) is the total number of degrees of freedom. Thus, the scaling 
factor for the velocities is 



A 



^ r \T(t)-l)' 



Within the Berendsen barostat the system is made to obey the equation of 
motion 

dp Po-P 

Tt = (5) 

where p is the value of the applied pressure, p is the instantaneous system 
pressure. In the isotropic case, the box volume is scaled by a factor 77, whereas 
the coordinates are scaled by 77 1 / 3 , where 

V = 1 ~ ^^(Po ~ p) (6) 

T 

and (3 is the isothermal compressibility of the system. The value of the coupling 
parameter r in our simulations was equal to 0.1 ps. 



3 Results 



3.1 Equation of State and Local Structural Properties 



A central place at the study of any phase transition or structural transfor- 
mations takes the consideration of the equation of states. Here, the equation 
of state was defined for the considered water model at the fixed temperature 
T = 290 K and different pressures. Results are shown in Fig. 2. As can be 
seen from the figure, the presented isotherm p{p) is well fitted by the cubic 
polynomial 

Pip) = c + c x p + c 2 p 2 + c 3 p 3 , (7) 
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Figure 2. (Color online) Main: Equation of state for liquid water at the temperature 
T = 290 K. Inset: Pressure dependence of the isothermal compressibility kt- 

Table 1 

Parameters of equation (7) of liquid water at the temperature T=290 K. 



Coefficients 


290 K 


c (g/cm 3 ) 


0.9996 


ci (g/cm 3 • Atm -1 ) 


6.204-10 -5 


C2 (g/cm 3 • Atm -2 ) 


-4.642-10 -9 


C3 (g/cm 3 • Atm -3 ) 


1.349-10 -13 



the values of the parameters Co, ci, C2, C3 are given in Table 1. The obtained 
equation allows one to find directly the isothermal compressibility kt as 




There is no any visible singularity in the density p(jp) as well as in the isother- 
mal compressibility kt{p) (see inset of Fig. 2) in the vicinity of the pressure 
p c ~ 2000 Atm, where LDL/HDL-structural transformations were observed 
[1]. This is direct evidence of the absence of the first-order-like transition at 
LDL/HDL structural transformations, while in amorphous ices this transition 
exists [27]. 

To perform a detailed study at microscopic level, we suggest to consider the 
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Figure 3. (Color online) Main: Pressure dependence of the local isothermal com- 
pressibility Kj, (left panel) and the difference (r 2 — ri) (right panel). Inset: Radius 
of the first two coordination shells vs. pressure. 

features upon the local spatial ranges around a water molecule within the 
boundaries of the closest coordination shells. For water the first two coordi- 
nation numbers defined as [38] 

j"T c 

N(r c ) = Aim / r 2 g(r)dr, (9) 
Jo 

take the values iV(ri) = 4 and N(r2) = 16, where r\ and r 2 is the radius 
of the first and the second coordination shell, respectively; n is the oxygen- 
number density. Therefore, on the basis of Eq. (9) the sizes of the first two 
coordination shells r\ and r 2 can be found for an arbitrary condensed phase as 
the distances, where N{r\) = 4 and N(r 2 ) = 16. Then, the local, short-ranged 
properties can be extended to these of the spherical layer enclosed between 
ri and r 2 of the volume Vq = (47r/3)(r2 — rf). The isothermal compressibility 
extended to the local volume will be defined as 



4 = ~TF ■ (10) 




Fig. 3 presents the pressure dependence of these characteristics: the local 
isothermal compressibility Kj,] the radius of the first two coordination shells 
ri and r 2 , and the difference (r 2 — ri) associated with the linear size of the 
considered local range. Contrary to the isothermal compressibility (inset 
of Fig. 2), which has a nonlinear monotonic decrease with pressure, the local 
isothermal compressibility demonstrates a significant jump at the pressure 
p c « 2000 Atm. A peculiarity in the vicinity of this value of pressure is also 
clear observed for the difference r 2 — r x presented in Fig. 3, although it is not 
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so obvious in the dependencies r\(p) and ^(p). Such result is a direct evi- 
dence that the structural transformations from LDL to HDL in water at the 
pressure p c ~ 2000 Atm are related with the structural rearrangement within 
the second coordination shell. Then, the next question arises naturally: Have 
these structural transformations an influence on the inherent microscopic dy- 
namics? One of the simplest way to answer on the question is to consider the 
vibrational density of states of the system. 



3.2 Vibrational Density of States 



The vibrational density of states (VDOS) characterizes the distribution over 
frequencies related with the molecular vibrational degree of freedom 



1 

2^ 



0(0) ■ 0(t) 



oo (tf(O) -0(0) 



dt 



(11) 



where 0(i) is the velocity of the mass center of a molecule at the time t. Fig. 
4 presents the defined VDOS of liquid water at the different pressures and the 
fixed temperature T = 290 K. As can be seen, VDOS has a bimodal form for 
all the cases. It is remarkable that the high-frequency peak position in VDOS 
is independent on pressure whereas the low-frequency peak changes both the 
intensity and the position u s . As can be see from the pressure dependence of 
ui s given in inset (a) of Fig. 4, the quantity u s has a non-monotonic step-like 
behavior in vicinity of the characteristic value of pressure p c ~ 2000 Atm. This 
indicates on the correlation of the low-frequency molecular vibrational proper- 
ties with the structural transformations at the pressure. While high-frequency 
mode corresponds to the inherent vibrational processes of the molecules, the 
time scale r s ~ l/ou s is related with the vibrational dynamics of the molec- 
ular conglomerates [40,41]. The immediate correlation of the low-frequency 
vibrations with the structural transformations at the local spatial scales in 
the range of the second coordination shell can be directly seen from inset 
(b) of Fig. 4. Recently, we have shown [1] that the increase of pressure gives 
rise the monotonic decrease of the tetrahedral order parameter, since the spa- 
tial structures of the five bonded water molecules changes. The low-frequency 
mode can be attributed to the vibrational motions of the separate molecu- 
lar tetrahedral complexes. So, one can concluded that the pressure-induced 
changes in the low-frequency mode is governed by the change of geometry of 
the local molecular conglomerates formed by five bonded molecules. 
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Figure 4. (Color online) Main: Vibrational density of states of liquid water at the 
temperature T = 290 K. Inset: (a) Pressure dependence of the low-frequency peak 
in VDOS indicated by arrow on main panel; (b) Radius of the second coordination 
shell T2 vs. characteristic time scale associated with the low-frequency vibrational 
dynamics, l/co s , for liquid water at the temperature T = 290 K. 

4 Summary 



Results of this study suggest explicitly the existence of the structural transfor- 
mations from low-density state to high density one at a "stable liquid region" 
of the water phase diagram. Namely, we found for water system at the temper- 
ature 290 K that the Amoeba model reveals the transformation in the vicinity 
of the characteristic value of pressure p c ~ 2000 Atm. It is necessary to note 
the corresponding structural transformation from LDL to HDL water is not 
marked in the temperature-pressure phase diagrams of water given recently, 
for example, in Ref. [2] (see Fig. 1) and Ref. [26] (Fig. 4). 

Strikingly, the transformation considered in this work has a lot of similar 
features with the polyamorphic phase transition from the low- to high-density 
amorphous forms of ice, which has been experimentally observed near 2000 
Atm and recognized as a first-order like transition (see work of O. Mishima 
et al. [27]). Furthermore, the transition in amorphous phase has a peculiarity 
to involve the movement of one water molecule from the second coordination 
shell to an interstitial site in the first coordination shell [22,42], that is very 
similar to LDL/HDL structural transformation considered here. Nevetheless, 
in the studied temperature-pressure range of stable liquid phase, there is no 
signatures of a phase transition in the thermodynamic sense between these 
two forms of water. 
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The isothermal compressibility kt of the system decreases steadily with pres- 
sure, while the local extension of the isothermal compressibility Kj. as well 
as other local structural and vibrational characteristics manifests a step-like 
behavior at the characteristic value of pressure p c m 2000 Atm. It is found 
that the LDL/HDL structural transformation in the vicinity of the pressure 
p c ps 2000 Atm has an influence on the low-frequency vibrational dynamics 
whereas the high-frequency molecular dynamics is independent on the pres- 
sure. Moreover, there is a strong correlation between low-frequency vibrational 
dynamics and the spatial scale of the second coordination shell. 
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